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SUMMARY 


A short theoretical study of aircraft aerodynamic model equations 
with unsteady effects is presented. The aerodynamic forces and moments 
are expressed in terms of indicial functions or internal state variables. The 
first representation leads to aircraft integro-differential equations of 
motion; the second preserves the state-space form of the model equations. 
The formulation of unsteady aerodynamics is applied in two examples. The 
first example deals with a one-degree-of-freedom harmonic motion about 
one of the aircraft body axes. In the second example, the equations for 
longitudinal short-period motion are developed. In these examples, only 
linear aerodynamic terms are considered. The indicial functions are 
postulated as simple exponentials and the internal state variables are 
governed by linear, time-invariant, first-order differential equations. It is 
shown that both approaches to the modeling of unsteady aerodynamics lead 
to identical models. In the case of aircraft longitudinal short-period 
motion, potential identifiability problems, if an estimation of aerodynamic 
parameters from flight data were to be attempted, are briefly mentioned . 
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INTRODUCTION 


One of the basic problems of flight dynamics is the formulation of 
aerodynamic forces and moments acting on an aircraft in arbitrary motion. 
For many years the aerodynamic functions were approximated by linear 
expressions leading to a concept of stability and control derivatives. The 
addition of nonlinear terms, expressing, for example, changes in stability 
derivatives with the angle of attack, extended the range of flight conditions 
to high-angle-of-attack regions and/or high- amplitude maneuvers. In both 
approaches, using either linear or nonlinear aerodynamics, it is assumed 
that the parameters appearing in polynomial or spline approximations are 
time invariant. However, this assumption was many times questioned 
based on studies of unsteady aerodynamics which go back to the twenties. 

A fundamental study of unsteady lift on an airfoil due to abrupt 
changes in the angle of attack was made by Wagner in reference 1. This 
work was extended by Theodorsen to computing forces and moments on an 
oscillating airfoil, whereas Kiissner and Sears studied the lift on an airfoil 
as it penetrates a sharp-edge or harmonically-varying gust, respectively 
(see reference 2). One of the first investigations of unsteady aerodynamic 
effects on aircraft motion was made by R. T. Jones in reference 3. He 
studied the effect of the wing wake on the lift of the horizontal tail. A more 
general formulation of linear unsteady aerodynamics in the aircraft 
longitudinal equations in terms of indicial functions was introduced by 
Tobak in reference 4. Later, in reference 5, Tobak and Schiff expressed the 
aerodynamic forces and moments as functionals of the state variables. 

This very general approach includes linear unsteady aerodynamics as a 
special case. A different approach to unsteady aerodynamics in aircraft 
equations of motion was introduced by Goman and his colleagues in 
reference 6. They used additional state variables, which they called 
internal state variables, in the functional relationships for the aerodynamic 
forces and moments. 

Despite the advancements of theoretical works, only a limited 
number of attempts were made to estimate aerodynamic parameters from 
experimental data and to demonstrate the importance of unsteady terms in 
aircraft equations of motion. In reference 7, a procedure for the estimation 


4 



of aerodynamic forces and moments from flight data was proposed. It 
starts with the estimation of stability and control derivatives. Then, the 
resulting residuals of the response variables are used for the estimation of 
unsteady terms. Reference 8 addresses identifiability problems for 
parameters in integro-differential equations. Examples of estimated 
indicial functions from simulated and flight data are given. Fourier 
functional analysis for unsteady aerodynamic modeling was applied to 
wind tunnel data of a triangular wing and a fighter aircraft in references 9 
and 10, respectively. It was shown that this modeling method was 
successful in computing the aerodynamic responses to large-amplitude 
harmonic and ramp-type motions. Finally, a concept of internal state 
variables for expressing unsteady aerodynamics was applied to wind 
tunnel oscillatory data and flight data in references 6 and 11. 

The purpose of this report is to summarize the approaches of 
references 5 and 6 to the formulation of aerodynamic model equations 
suitable for parameter estimation from experimental data. The report 
starts with expressing aerodynamic forces and moments in terms of 
indicial functions and internal state variables. Then, two examples of 
aerodynamic models for aircraft in small-amplitude motion are given. A 
discussion of these examples is completed by concluding remarks. 


AERODYNAMIC CHARACTERISTICS IN TERMS OF INDICIAL 

FUNCTIONS 

Using the results of reference 5, aircraft aerodynamic characteristics 
can be formulated as 

t , 

C a (t) = C a (0) + jc a( (t - t; i (t>) T ^ s ( r)dT (1) 

where 

C a (t ) is a coefficient of aerodynamic force or moment, 

£ is a vector of aircraft state and input variables upon which the 
coefficient C a depends, 

C a ( t ) is a vector of indicial functions whose elements are the responses 
in C a to unit steps in £ , and 

C a (0) is the value of the coefficient at initial steady-state conditions. 
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The indicial responses, C a , are functions of elapsed time (t - r) and are 

continuous single-valued functions of £{t). The indicial functions 
approach steady-state values with increasing values of the argument 
(t- r). To indicate this property, each indicial function can be expressed as 

C HJ (^ " ( T) ) = C HJ (°° ; ^ ( T) ) " F H,j ^ ~ V ^ { r) ) (2) 

where 

C a . (°°; £(r)) is the rate of change of the coefficient C a with £ j, in steady 

flow, evaluated at the instantaneous value of %j with the remaining 

variables £ fixed at the instantaneous values £(t) and 
the function F nr is called the deficiency function. This function 

approaches zero for (t - t) — » 

When equations (2) are substituted into equation (1), the terms involving the 
steady-state parameters can be integrated and equation (1) becomes 


c a (t) = C„(~; e «)) - JF 0£ (t - 1 ; £ tt)) T ^ £ ( rtdi 


(3) 


where 

C a (°°;£(£)) is the total aerodynamic coefficient that would correspond to 
steady flow with £ fixed at the instantaneous values £(t) , and 
F Q is a vector of deficiency fimctions 

If the indicial response C a ^ is only a function of elapsed time, equations (1) 


and (3) are simplified as 


c a (t) = C c (0) + J C Qi E (t - r) T £( T)dT 

= C a (0) + C B{ (~) T $ ( T) - 1 F a| a - t) t A 5 ( T)dr 


(4) 


When analytical forms of deficiency functions are specified, the 
aerodynamic model based on equations (3) or (4) can be used in the aircraft 
equations of motion for stability and control studies involving either linear 
or nonlinear aerodynamics. The resulting equations of motion will be 
represented by a set of integro-differential equations. 
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FORMULATION OF AERODYNAMIC FUNCTIONS USING INTERNAL 

STATE VARIABLES 


When the indicial functions are used in aircraft aerodynamic model 
equations, it is not clear, either from theory or experiment, what analytical 
form these functions should have. After postulating models for indicial or 
deficiency functions, questions about the physical meaning of terms in 
these models may still be asked. In order to avoid, at least partially, these 
questions, a concept of internal state variables for modeling of unsteady 
aerodynamics was proposed in reference 6. This approach retains the 
state-space formulation of aircraft dynamics, that is 

x = f(x(t),u(t)); ar(0) = jf 0 (5) 

by augmenting the aircraft states with the additional state variable i)(t). 
Then, the aerodynamic coefficients are formulated as 

C a (t) = C a (Z(t),ri(t)) (6) 

where 

^ ”£(*!(*), 5(f), 4(f)) (7) 

and 

5(*) = [*(f) T m( 0 t ] 

An example of equation (7) for a study of aircraft longitudinal 
dynamics is given in reference 6. Here, the internal state variable 
represents the vortex burst point location along the chord of a triangular 
wing. This location is described as 

7j*7+*?= no{a-T a a); |r/|<l (8) 

where 

t)q is the vortex burst point location under steady conditions, 

Tj is the time constant in the vortical flow development, and 
T a is the time lag in the same process caused by the angle-of-attack rate 
of change. 

The experimentally-obtained effect of the angle of attack and pitch rate on 
vortex point location is taken from reference 12 and is plotted in figure 1. 

The resulting curves were obtained by flow visualization on a delta wing 
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undergoing static test and forced pitching oscillations at two reduced 
frequencies. The effect of pitch rate is seen by comparing the dynamic 
vortex burst point location with the static point location. 


EXAMPLES 

The following two examples demonstrate the formulation of 
aerodynamic equations and equations of motion with unsteady 
aerodynamics. In these examples, only small-amplitude motion will be 
considered, thus leading to a system of linear equations. In the first 
example, aircraft one-degree-of-freedom (one d.o.f.) oscillatory motion about 
each of the three body axes is considered. The second example deals with 
short-period longitudinal motion. In both examples, the formulation of 
unsteady aerodynamics using indicial functions and internal state 
variables is considered. 


Harmonic Oscillatory Motion: 

In the development of aerodynamic models of an aircraft performing 
a one d.o.f. oscillatory motion, an approach using indicial functions and 
internal state variables will be considered. For the oscillatory motion in 
pitch, the functional relationships for the lift and pitching moments are 

C L (,t) = Ci(a(t),q(t)) 

C m (t) = C m (a(t),q(t)) 

Applying equation (4), the lift coefficient can be expressed as 

t j e t H 

C L (O = C i (0) + JC La (f-T)^-a(T)dr+^Jc Lg (t-T)^q(T)dT 
o T v o T 

£ J 

= C L (0) + C La (°o) a(t) - J F a (t - 1 ) — a( r)dr (9) 

0 T 

+ |rC L ,(«.) 9 (0-|:jF,((-r)^ :9 (T)dT 
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Similar expressions can be written for C m (t). Neglecting the effect of q(t) 
on the lift and taking into account only the increments with respect to 
steady conditions, equation (9) is simplified as 

C L (t) = CL a (o°)a(t)~ J F a (t - r)-^a(r)dr+^C Lq (oo) q (t) (10) 

For obtaining a model with a limited number of parameters, the indicial 
function is assumed to be in the form of a simple exponential 

C La (t) = a^l-e~ blt) j + c (11) 


Because 


lim Ct ( t) = a + c = Ci a (°°), 


equation (11) can also be written as 

CL a m = C La (~)-oe- b ‘‘ 

After substituting (11a) into (10) and applying the Laplace transform to 
equation (10), the expression for the lift coefficient is obtained as 


(11a) 


C L (s) = 


as t 


+ 77 C L -8 


a(s ) 


^ s + 6t V ^ 

\ X , 

where 

q(s ) was replaced by sa(s) and, for simplicity, Cl = Cl o (°°) and 
C L g =C Lq (oo). 

Using a complex expression for harmonic changes in a(t), that is 

a(t ) = oc^e 10 * = a^(cos((ot) + ism(cot)) , 
and replacing s by io), the steady-state solution to equation (12) is 


( 12 ) 


C L (t) = 


'L a 


a 




bj + co 2 


r lc a 6 , ^ 
■tr^La a ,2 , „2 


V b\ + (0 


a a sin(caf) 
a^acosicot) 


(13) 
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The introduction of reduced frequency 



and nondimensional time constant 


yields 

where 


*2 


_V^ 

bj£ 


C L ( t ) = Ci /a a a sin(ctf) + a^k cos (cut) 


Cl=C l (o o). 


■a 


fy 2 

l+rjk 2 




(oo)-a 


n 

1+T$k 2 


(14) 


(15) 


Similarly, the steady-state solution for the pitching-moment coefficient will 
he 


where 


C m (t) = C ma a^s\n{at) + C m ^a/Jicos{(Bt) 




(oo)-a 


rfo 2 

1 + t^k 2 



= C m (oo)- a 


l+r*k 2 


(16) 


(17) 


The parameters a and t 2 in equation (17) have, in general, different values 
from those in equation (15). 

When the internal state variable is used in formulating the unsteady 
aerodynamic effect, the development of a model for the lift coefficient starts 
with the equations 


C L {t) = Ci J (a(t),q(t),ri(t)) 

(18) 

T 1 f]+7j=T] 0 (a-T a d) 

(8) 
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For small perturbations, both equations can be linearized around a steady- 
state condition. Then, the linearized equations (18) and (8) will have the 
form 


C L (t) = C La a(t) + y C Lg q(t) + C h 77 (t) (19) 

T 1 f 1 +i 1 = -(T 1 +T cc )^6c (20) 

Applying the Laplace transform, these equations will be changed as 

C L (s) = C La a(s) + ^ C Lg q(s ) + C LjJ q(s) (21) 

(Tjs+Vjrjis) = ”( T i +r «)^ sa(s) (22) 

When equation (22) is substituted into (21) and q(s) is replaced by scc(s), 

C L (s) = C La cc(s) - ^ C Lti saM + |c V «(s) (23) 


Finally, introducing 

a = T l +T a&M c L and = TJ 1 
T x da ^ 

equation (23) will have the same form as equation (12). The preceding 
developments indicate that, for the indicial function given by equation (11a) 
and the internal variable given by equation (22), the model 

t , f 

Cl ( t ) = C La a(t ) - a J e~ bl{t ~ T) — a(r)dr +—C Lq q(t) 
o 

is equivalent to the model 

C L (t) = C La a(t) + ± C Lq q(t) + C Lv 


T 1 i]+r\ = -{T 1 +T a )^a 
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Models for an aircraft performing one d.o.f. oscillatory motion in roll 
and yaw can be developed in a similar way to that for the pitching 
oscillations. The rolling-moment coefficient is a function of the roll angle 
and rolling velocity 

Ci(t) = Ci{<l>(t),p(t )) (24) 

where the roll angle is related to the sideslip angle by the equation 

P = <f>sin(a) (25) 


For the indicial function 


< V ( > = C;/“>- ™ b,t 

the rolling-moment coefficient can be formulated as 

t , 

Ci(t) = C l (oo)p(t)-a\e- b ^ t - r) ^-p(r)dr+i-C l (~ ) p ( t ) (26) 

P J qt y P 

0 


which leads to its steady response 


where 


Q ( t ) = Cip <p A sin(a^) + Q <j) A k cos {(at) 


tV 

C ifi = C h (°°)sin(«) -a— sin(«) 

C; =Ci (oo)-a — ^-2 sin(a) 
p p 1+ifA: 2 


(27) 


(28) 


In the yawing oscillatory motion, the yawing-moment coefficient is a 
function of the yaw angle and its rate 

C n (t) = C n (yr(t),r(t)) (29) 

and the yaw angle is related to the sideslip angle as 

P = -y/cos(a) (30) 
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The yawing-moment equation takes the form 


C„(i) = C„/^«)-aje- 4 '«- l) ^j3(T)dr + lc nr (~)r«) 


( 31 ) 


and its steady response the form 


where 


C n it) = C n iff A sin(Grf) + C„ r \j/ A kcos(a)t) 


*2*2 

c np = c n p (°°)cos (a) - a - cos(oc) 


(32) 


(33) 


Qi r = Cn r (°°) + « - — ^72 cos(a) 

.L t r2 


For the interpretation of measured aerodynamic forces and moments 
in the forced-oscillation experiment, the model for an increment in the lift 
without any unsteady effect is usually postulated as (see reference 13) 

C L (t) = C La a(t) + -^;{c L(i a(t) + C Lq q(t)^ + {J^ C L ^q(t) (34) 

The unsteady version of the preceding equation will have to include two 
indicial functions, Ct (t) and Ct (t). Then the lift coefficient will be 

"a 

formulated as 


* J 

C L (t) = C L a(t)-\F a (t-r)-f-a(x)dr 
a i dr 


+ic,lU)-^F q (t-T)± g (T)dT 


(35) 


In both cases, the steady-state solution is given by equation (14) where, for 
the neglected unsteady aerodynamics, 


c La =c La -k 2 c L . 
Cl q =c Lq +c L . 


(36) 
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and, for the deficiency functions specified as 


F a = aje and F q = 056 

A 


-b<?t 


Ct ~ Ct — 

U a Lt a 




ai i + fy 2 a 2 i + ? 2 k 2 j 
r , r|fc 2 ^ 


( 37 ) 


a l 


i + ijk 2+a2 i + 4k 2 


From a comparison of equations (36) and (37), it can be concluded that the 

expressions in the parentheses are the unsteady counterparts to the 
derivatives Ct and Cj . . For large values of r and small values of k, the 

expressions in equations (37) can be simplified to those in equation (15). 
Similar comparisons can be made for the remaining aerodynamic 
coefficients. 


Short-Period Longitudinal Motion: 

The airplane short-period longitudinal motion can be described by the 
equations 

a = q+£®C z {a(t\q(t),S(t)) 

2m 

2 - (38) 

Zly 

In the following analysis, it will be assumed that the linear approximation 
to the aerodynamics contains only one unsteady term represented by the 
indicial function 


Cm a ~ ^m a (°°) ~ F a (t) 


(39) 
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Using simplified notation for the steady aerodynamic terms, the 
aerodynamic equations in (38) will have the form 


C z (t) = C Za a(t) + — C Zq q(t ) + C Zg 8(t ) 


d_ 

dr 

+ + C ms 8(t) 


C m ( t ) = C ma ait ) - j F a it - t ) — ai r)dr 


Specifically, for 


V <7 


Cm a it) = a(l~e blt) j + c 


the pitching-moment coefficient takes the form 

t , . 

C m it) = C ma ait)-a\e~ hl{t ~ x) — air)dr + — C mq qit) + C ms 5it) 


where 


= cait ) + a6j je blT ait - r)dr+ — C mq qit) + C ms 8it) 
0 

C m „ =a + c 


( 40 ) 


(41) 


Substituting (41) into (38) and introducing dimensional parameters, the 
equations of motion can be written as 

a = Z a a + Z q q + Z$8 

t 

q = Ca + B^e~ blX ait- r)dr + M q q + M§8 
0 

where the parameters in these equations are defined in table I. 
Introducing a new state variable 


(42) 


x a = je blX ait - r)dr 
0 

and the corresponding state equation for this variable 

x oc = a-b 1 x cc 
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equations (42) can be expressed in state-space form as 


a 


Nj 

o 
1 

a 


Z S 

q 

- 

C Mg B 

Q 

+ 

M 5 

x a _ 


, 

1 

o 
1 

x a _ 


0 


( 43 ) 


The characteristic polynomial of these equations has the form 

a=x 3 +k 2 i?+k 1 x+k 0 

where 

K 2 = -Z a - M q + bj 

Ki = Z a (M„ - b, ) - b : M q - CZ, (44) 

ir 0 =6 2 (Z a M,-Z g M a ) 


The state equations of the system under consideration can also be 
obtained by using the internal state variable defined by equation (20) as 


Tj+T, 


a 


77 = - 


= k a a-Ti l ri 




da 


(20a) 


When equation (20a) is combined with the equations of motion, the complete 
set of state equations is 


a 


Z a 

Z d 

0 

a 


■ 25 ' 

q 

= 

M a 

Mg 

Mg 

q 

+ 

M S 

V. 


_ka Z a 

k<x Z q 

-TJ 1 

_q_ 




(45) 


where M a and Mg are also explained in table I. After formulating the 

characteristic polynomial, it is found that its coefficients are equal to those 
defined by equations (44) for 


and 


a = 


Tl +Tg ^ d7? 0 
< Ti J da 



bj = TJ 1 
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It can, therefore, be concluded that equation (43) and (45) represent the 
same dynamical system. As in the previous example, the system 
description using either indicial functions or internal state variables can be 
identical for specific forms of indicial functions and equations for internal 
state variables. 

The preceding development shows that the introduction of one 
indicial function of the form specified by equation (39) into the aerodynamic 
model equations results in the increase of the order of the characteristic 
polynomial from two (no unsteady aerodynamics) to three. Any further 
addition of indicial functions into equation (42) means an additional 
increase in the order of the characteristic polynomial by one. From a 
simple observation of equation (43) or (45), it is also evident that it is not 
possible to estimate all the parameters in these equations from the 
measurements of a(t), q(t), and 5(t). To assure parameter identifiability, 
equation (43) would have to be transformed into a canonical form proposed 
in reference 14. 

In stability and control analysis where no unsteady aerodynamics is 
considered, the pitching-moment coefficient is formulated as 

C m =C ma a + ^ (C m . a + C m<? g) + C mg S 


It is expected, therefore, that the integral in equation (40) 

t 


I = \F a {t-x) — or(T)dr 
i dr 

0 


(46) 


f £a) 

should be a counterpart of the term I — \C ma • Th e reduction of this 

integral to the a -term can be demonstrated by approximating a(t) by a 
Fourier series 


oc(t) = Aq + (Aj - iBj )e io * + (A 2 — il $2 +... 

which leads to 

a(t ) = i(o{Ai -iBi)e iwt +i2o)(A2 -iB 2 )e i2at +... (47) 
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Substituting (47) into (46) results in 


t 

I = ico{A 1 -iBi)e ia * j F a (r)e- io}T dT 

0 

t 

+ i2(o{A 2 -iB 2 )e i2cat \ F a ( T)e- i2o)r dT+... 

0 


( 48 ) 


The exponential functions in (48) can be further expanded in exponential 
series 


e im _ i + i m + -(i (0 ' U ) 2 +" W 

2 

e l2m _ i + i2m+2(i(or) 2 +... 


In order to maintain the approximation of the integral to the first order in 
frequency, it is sufficient to consider only the first terms in the exponential 
series. Then, all the integrals in (48) will be the same and equation (46) can 
be simplified as 

t 

I = a\F a { r)dr (49) 

o 

As a result of this simplification, the counterpart of is proportional to 

the area of the deficiency function. A similar conclusion is stated in 
reference 4 for simple harmonic motion of an aircraft. 

For a demonstration of aircraft longitudinal motion with and without 
unsteady aerodynamic terms, equations (42) and their simplified version 

cc = ZfyOC + Z a q + Zx5 

9 (50) 

q = M a a + M^a + M q q + M$S 

were used. Aircraft characteristics and flight conditions are summarized 
in table II. The unsteady parameter bj was selected as bj = l|sec -1 j which 

corresponds to the nondimensional time constant Tj =51.3. The parameter 
a was evaluated from the relationship between the derivative C m& and the 

area of the deficiency function 
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as a = 0.05. Because C ma = a + c , the parameter c = -0.23. 

In table III, the computed damping coefficients and frequencies of 

motion from equations (42) and (50) are presented. The values of these 
parameters indicate that the replacement of the terms C mQ a and C m ^ a by 

the indicial function C ma ( t ) has a negligible effect on the damping 

coefficient and only a small effect on the frequency. Figure 2 shows the 
computed time histories a(t) and q(t) for the given input 8{t). As could be 
expected from the results in table III, the output variables for both cases 
differ only slightly. Small differences in a{t) and q{t) might indicate 
possible problems when estimation of unsteady parameters from flight data 
is attempted. 


CONCLUDING REMARKS 

A short theoretical study of aircraft aerodynamic model equations 
with unsteady effects is presented. First, the aerodynamic forces and 
moments are expressed in terms of indicial functions. This formulation 
can be modified by including steady values of aerodynamic coefficients, 
corresponding to instantaneous values of state and input variables, and the 
so-called deficiency functions. A deficiency function defines the difference 
between the indicial function and its steady value. When the concept of 
indicial or deficiency functions is used, the resulting aircraft model is 
represented by a set of integro-differential equations. In the second 
approach to the modeling of unsteady aerodynamics, the so-called internal 
state variables were used. These variables are additional states upon which 
the aerodynamic coefficient depends. Modeling based on internal state 
variables preserves the state-space representation of the aircraft equations 
of motion. 

The formulation of unsteady aerodynamics is applied in two 
examples. In these examples, only linear aerodynamics are considered 
thus limiting the application to aircraft small-amplitude motion around 
trim conditions. In order to further simplify the aerodynamic model 


19 



equations, the indicial functions are postulated in a simple exponential 
form and the internal state variables are governed by linear, time- 
invariant, first-order differential equations. 

In the first example, a one-degree-of-freedom harmonic motion about 
one of the aircraft body axes is considered. In the second example, a 
longitudinal short-period motion is studied. In both examples, it is shown 
that the formulation using either indicial functions or internal state 
variables leads to identical models. Further, it is shown that the unsteady 
terms in the models are the unsteady counterparts of the aerodynamic 
acceleration derivatives. From an observation of the developed longitudinal 
equations of motion, it is evident that it will be impossible to estimate all 
aerodynamic parameters from measured input/output data. In addition, a 
simple numerical example of the short-period motion of a fighter aircraft 
indicates only small differences in the output time histories with the 
unsteady effects being either included or ignored. These small differences 
might create further problems when estimation of unsteady parameters 
from flight data is attempted. 
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Table I. - Definition of parameters in equations (38) and (39). 


Z a = psv c z 
a 2m 

A/f pV 2 Sc 

Mry = C m 

2/y “ 

Z * = 1+ 4 S >< 

M a = pVSl ' C m 
4/y a 

7 -P sv r 

Zs - 2m ° Z ‘ 

M q - pVS? C m 

q 4/y 

C = pV2S 'c 

2/y 

M„ = pV2S ~ C C m 

1 2/y m 1 

B = pV2S 'a bl 

2/y 1 

M pV 2 Sc _ 

M<5_ 2/y Cm « 


Table II. - Characteristics of an advanced fighter aircraft 
and flight conditions. 


c = 3.51 m 
<S = 37. 16 m 2 
m = 15000 kg 
I Y = 170000 kg-m 2 
p = 0.56 kg/m 3 
V = 90 m / sec 


C z =-2.7 

C z =-36. 

<] 

C Zs =-0.83 
C m „=-0.18 
= —2.5 
Cm q =~ 10 . 

C ms =-0-88 
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Table III. - Damping coefficients and frequencies from simulations with 

and without unsteady effects. 



damping coefficient 

frequency 

with unsteady effects 

0.4859 

0.6317 

without unsteady 
effects 

0.4979 

0.5953 



Figure 1. - Variation of internal state variable with angle of attack 
in static and oscillatory tests. 
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